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By means of the continuous unitary transformation similar to a general scheme of the Renormal- 
ization Group (RG) procedure we study the issue of symmetry breaking and pairing instability in the 
system of interacting fermions. Constructing a generalized version of the Bogoliubov transformation 
we show that formation of the fermion pairs and their superconductivity/superfluidity can appear 
at different temperatures. It is shown that strong quantum fluctuations can destroy the long-range 
order without breaking the fermion pairs which may still exist as incoherent and/or damped entities. 
Such unusual phase is characterized by a partial suppression of the density of states near the Fermi 
energy and by residual collective features like the sound-wave mode in the fermion pair spectrum. 



Formation of the fermion pairs is a common phe- 
nomenon for various physical systems of the interact- 
ing particles such as: electrons, nucleons, atoms and 
quarks. Binding energy and the spatial extent of fermion 
pairs may vary from case to case depending on particular 
species and on specific interaction mechanism. To give 
some examples let us mention that pairing can be driven 
by: 

(i) exchange of phonons (in classical superconductors, 
MgB 2 , etc), 

(ii) exchange of magnons (superconductivity of the 
heavy fermion compounds), 

(hi) strong correlations (the high T c superconductors), 

(iv) Feshbach resonance (superfluidity of the ultracold 
fermion atoms), 

(v) or by other effects (nucleon pairing in nuclei, su- 
perfuidity of the neutron stars). 

Usually formation of the fermion pairs goes hand in 
hand with appearance of the order parameter which con- 
sequently leads either to superconductivity (for charged 
particles such the conduction band electrons or holes) or 
to superfluidity (for electrically neutral objects like 3 He 
or the ultracold fermion atoms in magnetooptical traps) . 
However, a simultaneous formation of pairs and emer- 
gence of the symmetry broken phases needs not be a rule. 
We will show here example that both these phenomena 
are distinct and happen to coincide at the same critical 
temperature T c only when the quantum fluctuations are 
weak. 



I. HAMILTONIAN OF THE INTERACTING 
FERMIONS 

System of the interacting fermions can be described by 
the following Hamiltonian 
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where is a single particle energy for a given momentum 
k and a corresponds to additional quantum numbers like 
for instance spin f, j for electrons, the angular momen- 
tum for atoms or the isospin for nucleons. The two-body 
interactions are described by the second term with the 
potential E/k,k'(q)- We use in the standard notation 
for the creation (annihilation) operators c k (cko-)- 

In general there can arise various kinds of ordering, for 
instance: ferromagnetism, antiferromagnetism, charge 
ordering, superconducting BCS state, etc. We will fo- 
cus here on the pairing instabilities. For this purpose we 
further consider the Hamiltonian reduced only to q = 
channel 
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with £k = £k~ M an d we assume the two-body potential to 
be attractive Vk,k' < 0. We will investigate this reduced 
BCS Hamiltonian using the nonperturbative method (de- 
scribed in section III) which belong to a family of the 
Renormalization Group techniques 



II. RENORMALIZATION GROUP APPROACH 

Thermodynamics of the system (total energy, specific 
heat, pressure, etc) can be computed from the partition 
function Z — Tre~ H / kBT . It is convenient to express 
Z using the path integrals over the Grassmann variables 
ipk, a, "0k o- (which formally represent the eigenvalues of 
the annihilation Ck,<r and creation c k operators) 



Z = J D[ip,if>*] e~ s . 
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The action consists of two contributions S = So + S7 
where the quadratic term 



So = E / (V„-£k) lpk,<r 

~ J k 



(4) 



corresponds to a free part and integration in the equa- 
tion Q runs over the four- vector k = (iw n ,k) with the 
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fermion Matsubara frequencies u) n = (2n+ l)7rfcsT. The 
second quartic term refers to the two-body interactions 
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As far as the dynamic quantities are concerned (various 
correlation functions) they are derivable directly from the 
generating functional 
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where Xk,a and xt a are the Grassman source 
fields. For instance, the single particle excitations 
can be determined via the two-point Green's function 

«Xk,„ SXk.v G[X>X L=o,x*=o- 

Physical properties of the system under consideration 
depend predominantly on such fermion states which are 
located near the Fermi surface. Other states distant from 
the Fermi energy are less relevant therefore it is useful 
to make a distinction between their contributions to the 
partition function 



£> <A [^,V>*] / D >A [ip,<ip*] eW'W 
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where symbol D <A [tp,ip*] corresponds to the fermion 
states whose distance from the Fermi energy is smaller 
than a given cut-off \ek — £k F | < A. In the numerical 
RG method pj one first integrates out the high energy 
fermion states. After completing such integration one is 
left with only the low energy states. Partition function is 
then Z = J D <A [^,^*} e~ s ^^ where the renormal- 
ized action is defined by 



(8) 



This action JHJ) can be then cast into © simplifying the 
calculation of the generating functional 



Q[x,x*} = 



log 



(9) 
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Such idea of mode elimination has been introduced in the 
theoretical physics a long time ago @- Wilson adopted 
it to the solid state physics by proposing a sequential 
integration of the fermion fields down to some small cut- 
off A. Reducing bit by bit the cut-off A to infinitesimally 
small values he was able to study emergence of the critical 
phenomena For a more specific discussion of various 
RG formulations we recommend the the review papers 

a 

However, in the case of symmetry broken phases (such 
as the superconducting state) the simple scaling proce- 
dure usually fails due to a natural lower boundary cut-off 



(the energy gap A of the single particle excitations). It 
has been even claimed that conventional RG techniques 
are blind with respect to the symmetry-broken states Q . 
The situation is not that much hopeless, there are pos- 
sible routes to circumvent this problem. Let us mention 
two of them 

(i) one can impose by hand an infinitesimal symmetry- 
breaking parameter A (A) at a certain initial con- 
dition A = Ao and then its physical meaning would 
eventually establish upon the flow A — > 5] , 

(ii) in order to eliminate the interacting part Si one 
can enlarge the Hilbert space by introducing the 
linear coupling to the boson fields <E>, <&* via the 
Hubbard - Stratonovich transformation 0. 

Using the first procedure one must be cautious how to 
constrain the small symmetry breaking parameter while 
in the latter method there rise additional complications 
dealing with the new boson fields. Perhaps the second 
option is more natural because after all the interactions 
are always mediated by certain boson fields (phonons, 
photons, gravitons, etc). Unfortunately in practice it is 
hard to go beyond a simple saddle point approximation 
except than by taking into account the small Gaussian 
corrections around it. In the next section we shall briefly 
describe a more convenient procedure which has been 
proposed in 1994. 



III. CONTINUOUS CANONICAL 
TRANSFORMATION 

An alternative approach to deal with the many-body 
effects (which in particular is suitable for studying the 
symmetry broken phases) has been invented by Weg- 
ner pfl and independently by Wilson and Glazek [|J. 
Instead of integrating out the high energy states (fast 
modes) one constructs a continuous unitary transforma- 
tion H(l) = U{l)HW (I) with a purpose to simplify the 
Hamiltonian either to diagonal or at least to a block- 
diagonal form. This is achieved through a series of in- 
finitesimal steps H — ► H(l) — ► H(oo), where I stands 
for some formal continuous parameter. In a course of 
transformation the Hamiltonian evolves according to the 
differential flow equation 



8 l H{l)= fj(l),H(l) 



(10) 



where the generating operator is defined as fj(l) — 
-U(l) di&(l). 

This method has a similarity to the traditional RG 
scaling procedure because 

(i) diagonalization of the high energy sector occurs 
mainly during initial part 

(ii) while the low energy states are worked out at the 
very end of transformation. 
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Roughly speaking, a role of the Wilson's cut-off energy 
A is played here by A' = 

Advantage of such new procedure becomes particularly 
clear when investigating the mutual relations between 
the high and small energy states (i.e. between the fast 
and small modes). In the continuous canonical transfor- 
mation one treats both energy sectors on equal footing 
throughout the whole transformation process. Thereof 
their feedback effects are feasible to study. 

For the Hamiltonians H — H + Hi (where H is a 
diagonal part in a given representation) it has been pro- 
posed to choose 



rK0= #o(0,#i(0 
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so that the off-diagonal term is eliminated at the asymp- 
totic point 



lim Hx{l) = 0. 
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There are also other possible ways for constructing the 
generating operator fj. Some a survey we recommend the 
recent review papers 0, El . 

As far as the diagonalization is concerned the continu- 
ous canonical transformation turns out to be a rather 
convenient tool. However, if one needs the correla- 
tion functions then the situation becomes more cumber- 
some. For example, to determine the correlation func- 
tions (A(t)B(t')) one needs the statistical average 



Tr 



= Tv{u{l)e-^6U\1)) 



Tr 
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where 0(1) = U(V)OW(l). The easiest way to carry out 
the statistical averaging l|13|) is in the limit I — > oo 
when ff (oo) becomes (block-)diagonal. However, this re- 
quires a simultaneous transformation of the observables 
O — > 0(1) — ► O(oo). The corresponding flow equation 
is given in a familiar form 



diO{l)= f)(l),0(l) 



(14) 



Thus, calculation of the correlation functions is rather 
more relative to the projection techniques. 

IV. THE BILINEAR HAMILTONIAN 

To illustrate how the flow equation method actually 
works we first consider the bilinear Hamiltonian 



k 



which can be solved exactly, in particular by a single 
step Bogoliubov transformation 11]. The off-diagonal 
terms in the Hamiltonian (|15fl can be thought as resulting 



from the mean field approximation for the weak pairing 
potential 14c, k' < with a usual definition of the order 
parameter A k = -J2k' ^k,k'(c_k'j.Ck'T)- 

To get the rigorous solution we now construct a con- 
tinuous transformation such that the convoluted states 
|k, t> and | — k, J.> will be decoupled. This continuous 
process depends on a distance from the Fermi surface 
|£k| (and the same holds for computation of the coher- 
ence factors tik, «k)' 

Using the Wegner's proposal (|llf> we obtain the two 
coupled flow equations 



di &(/) 
di A k (Z) 



4£ k (0|A k (0| 2 
-4(£f k (0) 2 A k (0 
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for /-dependent quantities £k(0 and A k (Z). The second 
equation l|17|) yields 



|A k (0| = |A k |exp^-4 / dl'[^(l')} 



(18) 



which proves that indeed the off-diagonal terms gradually 
diminish under the flow I — ► oo. There is a singular point 
k = kp which is unaffected by the transformation but in 
the thermodynamic limit (i.e. for a macroscopic number 
of particles N) its role becomes marginal. 

By combining the equations (|16I17|1 one gets the fol- 
lowing invariance 



5 z {(c; k (0) 2 + |A k (0| 2 }=0. 
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which implies that in the limit / — * oo the eigenvalues 
have a character of the Bogoliubov spectrum fl3| 



£ k (oo) = sgn(a)V(a) 2 + |A k | 2 . 



(20) 



In figure ^ we illustrate evolution of the parameter 
A k (Z) which initially was assumed to be constant A k = 
A. In the first turn disappearance of A k (Z) occurs for 
states distant from the Fermi energy and then in the sec- 
ond turn to states located nearby to fi. This evolution is 
accompanied by the renormalization of fermion energies 
£ k (0- Finally for / — > oo the quasiparticle energies evolve 
to the gaped Bogoliubov dispersion (J2DJ. 

In order to specify the complete single parti- 
cle spectrum we now derive the correlation function 
(c ko -(i)c k(T (i')) where time evolution is given by the stan- 



dard relation 0(t) 
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As emphasized in 



the previous section this requires determination of the 
Z-dependent single particle operators. From a detailed 
analysis [l2T ] we obtain the Bogoliubov-type parameteri- 
zations 
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FIG. 1: Variation of A^(l) (panel on the right) and £k(0 
(panel on the left) versus a distance from the Fermi surface 
£k F = 0. We assumed a constant parameter Ak(0) = A 
in the model Hamiltonian l|15^ . Energies presented in this 
figure are expressed in units of A and the flow parameter in 
units of A -2 . Note that for a given I the renormalizations are 
practically completed outside the energy window of the width 
A 1 — as marked by the arrows. 



with the initial values Uk(0) = 1, «k(0) = 0. From itT^ 
we derive the following differential equations for the co- 
efficients Uk(0; u k(0 

diu k (l) = 2£ k (Z)A k (/) v k (l), (23) 
d l wk(i) = -2£ k (/)A u (/) u k (l). (24) 

It can be easily checked that equations (|2^1l2^|) lead to the 
following invariance of Z-dependent coefficients |itk(Z)| 2 + 
bk(0| 2 = 1- This invariance assures that the fermion 
anticommutation relations are obeyed for arbitrary level 

of the transformation |ck CT (0j CicV'(0} = ^M'^cr.a'- 

We have previously shown ^| that the asymptotic 
I = oo values of the coefficients Uk(l) and Vk(l) coincide 
with the usual BCS factors 



«k(oo) 

M k (oo)w k (oo) 



2^ 



(25) 
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and v^(oo) = 1 — u^(oo). In figure |21 we show the l- 
dependent factors u^(l), v^(l) versus the energy mea- 
sured from the Fermi surface. We can notice that for ar- 




FIG. 2: The BCS coherence factors u k (Z) and u k (Z) versus the 
energy £ k for some representative values of Z as indicated in 
the legend. 



bitrary Z the asymptotic values are reached by all fermion 
states located outside the energy cut-off A' = 

The complete single particle excitation spectrum is 
made of the particle and hole contributions. The cor- 
responding spectral function is found |l2j | to be 

A(k, w) = u£(oo)<S(w - fk(oo)) + vl(oo)S(ui + &(oo)) 

(27) 

In a straightforward manner we can arrive at the follow- 
ing expressions for the average occupancy of k momen- 
tum and for the order parameter 



(CkT^k?) 
(c-kiCkt) 



1 - 



6 



-tanh 



£k (oo) 



£k(oo)~2fc B T 
A k , , Ck(oo) 



26c (c 



-tanh 



2k B T 



(28) 
(29) 



The Bogoliubov spectrum (|20|l together with the expec- 
tation values expressed in equations I|28I29|I reproduce 
the rigorous solution for the bilinear Hamiltonian 1)150. 



V. PAIRING IN THE STRONGLY 
CORRELATED FERMION SYSTEM 

The two particle interactions in JQ) or in the reduced 
Hamiltonian J2J can be canceled out exactly by introduc- 
ing the Hubbard-Stratonovich fields $, $*. Skipping a 
detailed derivation we assign here the following effective 
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Hamiltonian to the resulting fermion and boson degrees 
of freedom Q 



-7= E [ fe qC-kaCk+q, T + k' 

k,q 



(30) 



In the simplest step one can study ((3U(I on the level of sad- 
dle point approximation when boson fields are replaced 
by constant number. This is a mean field approxima- 
tion and formally it is equivalent to linearization of the 
boson-fermion term via 

& q C_ k)i C k+q , T ~ (S^)c_ kli C k+q , T + S^(c_ k)i Ck+q,T) 

- (S q ) (c_ kii c k+q , T ). (31) 

When substituting 131( 1 to the Hamiltonian (|30|) the 
fermion and boson subsystems become decoupled from 
one another. The physical fermion part acquires then 
the bilinear structure 1 (15(1 with the parameter A k = 
^k,o(&q=o)- The mean field properties of the model l|3T)|l 
have been summarized in the review paper |l5j . 

In this section we show how to go beyond the mean 
field scenario using the flow equation method. Our main 
strategy is to design a continuous canonical transforma- 
tion U(l) which, step by step, dismantles the fermion 
from boson degrees of freedom. Technical remarks con- 
cerning such transformation has been given in our pre- 
vious work where we constrained the Hamiltonian 
into H — H + Hb-f with H — Hp + H B denot- 
ing the independent fermion and boson contributions 
and Hb-f corresponding to their interaction. We have 
followed the idea proposed by Wegner Q by choosing 
77(0 = [Hq{1),H B -f(1)} which yields 



^) = -^E[^)- e -kW-£k +q (0] 

k,q 

x fu k)q (J)&qC_ k j.e k +kT - h.c.J • (32) 



All the Z-dependent quantities have been determined 
by us selfconsistcntly using the iterative numerical Rungc 
Kutta algorithm to solve the set of differential flow equa- 
tions (16-21) presented in the Ref. We studied a 
situation with the fixed number of particles in the system 
ntot = E k , CT (cLck<r) + 2£ (S q S q ). Among the impor- 
tant physical effects obtained for the fixed point I — > oo 
we could point out that: 

(c) boson particles (which can be thought as the 
fermion pairs) acquire a finite mobility owing to 
the interaction with the itinerant fermions, 

(b) fermions in turn are affected by the boson particles 
and this effect shows up by a loss of the single parti- 
cle states near the Fermi surface (such depletion of 
the density of states is often called in the literature 
as pseudogap), 



(c) in addition to the pseudogap there appears a reso- 
nant scattering between fermions (17|. 

Let us remark that the resonant-type scattering processes 
have been known for a long time in the nuclear physics. 
Such Feshbach mechanism is nowadays widely explored 
experimentally and theoretically for the atomic gasses 
cooled to ultralow temperatures enabling observation of 
such quantum phenomena like the Bosc Einstein conden- 
sation (BEC) and ultimately also the atomic superfuidity. 

Since fermion states located near the Fermi energy get 
combined with the boson species it is natural to expect 
that the single and two-particle properties are going to 
affect each other. We studied systematically their in- 
terplay within the flow equation method. To derive the 
single particle spectrum we had to transform the annihi- 
lation c k(T and creation c k(T operators using ((32(1 . From 
the flow equation ((14(1 we derived the following general- 
ized Bogoliubov transformation fl8| 



Mi 1 ) 

^ ki (0 



— y 



q#0 



Wk(0 

Pk,q(0 
t,q(0 



Ck| 



Mi) 



--kj. 



q C q+kT 



r k,q(0 



V-q-ki 



(33) 



In consequence the single particle spectral function was 
found to have a different structure than the usual 
BCS result l|27|l. For temperatures T < T c it con- 
sists of two narrow quasiparticle peaks at energies u) — 
i\/ (^k — m) 2 + an d a certain amount of the damped 
states forming an incoherent background outside the su- 
perconducting gap (see the left panel in figure OJ). When 
traversing the critical temperature T c to a normal state 
the gaped Bogoliubov-type spectrum seems to be pre- 
served, however for increasing temperature the shadow 
branch becomes more and more damped (see the right 
h.s. panel in figure 0). Physically it means that fermion 
pairs have no longer an infinite life-time. Finally, for tem- 
peratures exceeding a certain characteristic value T* the 
Bogoliubov modes are completely gone and there remains 
only a single peak at the renormalized energy £ k (oo) |18| . 

To check a direct impact of the above mentioned be- 
havior on the pair correlations we investigated the fol- 
lowing correlation function 

^(c_ ki (t)c q+kT (i) c q+k , T (i')cl kU (t'))- (34) 

k ,k' 

From the flow equation 11411 for the pair operators 
Sk c-kj.c q +kT we obtained the corresponding spectral 
function 

A vair (q, u>) = A/q 8 (u - £ q ) + A^ c (w) . (35) 

It contains the quasiparticle peak at w = Eq and the 
incoherent background A™ c (u>). For T < T c the quasi- 
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FIG. 4: Derivative of the fermion pair dispersion dEq/dq x 
(where _E q = _E q (oo)) for the normal (the dashed line) and 
for the pseudogap state close to T c (the solid curve). Note a 
tendency for a qualitative changeover from the parabolic to 
linear relation of _E q (oo). 




FIG. 3: The effective single particle spectrum of fermions 
A F (k, uj) for the superconducting state (panel on the left h.s.) 
and in the normal state near the critical temperature T > T c 
(panel on the right h.s.). 



particle peak is well separated from the incoherent back- 
ground and, in the limit q —> 0, has the important sound- 
wave dispersion _E q = c |q| (like the collective sound-wave 
branch in the superfhiid state of 4 He). Unfortunately in 
the case of charged fermions such as the conduction band 
electrons this mode is usually pushed to the huge plasma 
frequency because of the strong Coulomb repulsion. For 
the electrically neutral objects this Goldstone mode is a 
hallmark of the symmetry broken phase. 

Above the transition temperature and close to T c there 
are still some residual collective features possible to ob- 
serve. At small momenta the qusiparticle peak lu = 
overlaps with the incoherent background therefore the 
long wavelength limit is not suitable for appearance of 
the collective effects (strictly speaking for q — > the 
Goldstone mode is replaced by a parabolic dispersion). 
A remnant of the Goldstone mode splits off from the in- 
coherent background for finite momenta above a certain 
critical value q crj ;t- This is illustrated in figure 01 

Collective features seen in the pair-pair correlations (or 
in the density-density correlations) tell us directly about 
a presence or absence of the long range order which is 
necessary for the onset superconductivity/superfluidity. 
As a matter of fact the single particle properties are 



then ill defined (at least in a vicinity of the Fermi en- 
ergy). Instead of single fermions one should rather 
think about fermion pairs as good quasiparticles. Al- 
though the fermion pairs are necessary for superconud- 
ctivity/superfluidity the other way around this is not 
valid. It has been shown by Eagles ,19], Nozieres and 
Schmitt-Rink 20] and by a number of other authors that 
value of the binding energy (magnitude of the gap in the 
single particle spectrum) does not scale linearly with the 
superfluid phase stiffness n s which determines the tran- 
sition temperature T c . Our work indeed confirms that 
existence of the preformed (i.e. incoherent) fermion pairs 
is a natural expectation while approaching T c from above. 
Such type of situation can be encountered in the high T c 
cuprate superconductors where the strong quantum fluc- 
tuations are driven by the unscreened Coulomb repulsion 
between electrons and by the reduced dimensionality of 
CuC-2 planes |2l| . Pseudogap have been also unambigu- 
ously observed in the ultracold fermion atoms close to 
the unitarity limit (i.e. on the Feshbach resonance) [22|. 



VI. SUMMARY 

We have presented the method of continuous unitary 
transformation originating from a general scheme of the 
RG scaling procedure. This non-perturbative technique 
overcomes usual problems of the standard RG methods 
in application to the symmetry broken states. We have 
illustrated it on the exactly solvable case of the bilinear 
Hamiltonian. This new method is moreover capable to 
study the possible feedback effects between the fast and 
slow modes treating both sectors simultaneously through- 
out the whole continuous transformation. 

Applying this method to the strongly interacting 
fermion system we have shown that formation of the 
fermion pairs need not be accompanied by the transi- 
tion to superfluid/superconducting state. Strong quan- 
tum fluctuations can suppress the long-range coherence 
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(ordering) so that effectively the fermion pairs exist even 
in the normal state above T c . Such preformed fermion 
pairs could be observed experimentally by e.g. probing 
the single particle spectrum (using the STM or ARPES 
spectroscopies) or in measurements of the correlations 
between pairs (via any experimental technique sensitive 
to the pair susceptibility). 
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